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Abstract 

We find, within the pQCD dipole model for DIS processes, a rapid increase with energy of the 
scale of parton momenta in the current fragmentation region in the limit s — > cxo and fixed . 
We explain the equivalence between the dispersion representation over for LT pQCD zero angle 
amplitude of 7* + T ^ 7* + T scattering at large energies and the kt factorization, and use it to 
evaluate the scales of the hard processes in the current fragmentation region. We derive within 
the black disc (BD) regime the modified Gribov formula for the total cross-section of the DIS. We 
evaluate the coherence length of the processes relevant for the BD regime and find that it increases 
with energy as s^'^ i.e. significantly slower than in the parton model {l/2mNX - the loffe length) 
as well as in pQCD. In the BD regime we estimate the gluon densities, local in the coordinate 
space, and find that they do not decrease with the energy. We discuss briefly how the new pQCD 
phenomena may reveal itself in the proton-proton collisions at LHC. 



* Talk given by B.Blok at the small x physics 2008 conference at Columpari, Crete, July 2008. 
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I. INTRODUCTION 



The aim of this talk is to demonstrate the existence of the new hard QCD phenomena 
in DIS at small x = Q'^/s and large but fixed Q^, in the vicinity of the onset of the BD 
regime, in the kinematic domain where the pQCD is still valid, and to apply them to the 
physics to be observed at the LHC. We use the dipole model and kt factorization 0] as 
a tool and find that in the limit xb —>■ and large but fixed the essential transverse 
momenta of partons within the dipole are rapidly increasing with energy in contrast with 
the limit Q"^ ^ oo and fixed xb- In the case of the longitudinal photon cross-section ai 
the essential region of integration over the transverse momenta of the constituents within 
the quark dipole is k'f ^ Q"^/^ at moderate and large x but increases with the energy at 
sufficiently small x. Similar effects occur in the case of the transverse photon cross-section 
(Tt, where the relevant invariant masses increase with the energy even more rapidly. This 
result is in a sharp contrast with the regime of moderate energies like the energies of HERA 
(s < 10^ GeV^), where dominant configurations in the cross section qq have invariant masses 
M\qq) < Q\ 

The rapid increase of the characteristic transverse scales in the fragmentation region has 
been discussed in ref. jsj, however within the black disk (BD) regime. Our new result is the 
prediction of the increase of the jet momenta in the fragmentation region, in the kinematical 
domain where methods of pQCD are still applicable. To visualize analytically the origin of 
this phenomenon in pQCD we use the double logarithmic approximation which is a rather 
rough approximation. The calculations show that in ai the relevant invariant masses of the 
q — q pairs (dipoles) increase to ~ 1.5 — for the quark dipole- nucleon interactions 
at invariant energies s = 10^ - lO^GeV^. Let us stress that this new regime is characterized 
by the increase of the parton transverse momenta with energy in the photon fragmentation 
region and thus it is different from the well known phenomena of the diffusion to the large 
transverse momenta in In plane in the center of rapidity . 

We have carried our calculations in the leading logarithmic approximation. However, we 
expect that qualitatively our results will be valid in the NLO and beyond, because to a 
large extent the NLO corrections are taken into account in the phenomenological fits that 
fix the parameters of the parton distributions in the LO approximation. Remember that 
the parton distributions have been extracted from data in the limit of — oo and fixed x 
where effects discussed in the paper should be a correction. 

For the deep inelastic lepton scattering off a nucleon we evaluate the dependence on the 
energy s and transverse momenta of the three characteristic scales. The first scale is 
the maximum in the distribution of jets in the current fragmentation region of the total 
cross-section over -M^ . Then we evaluate the region of that gives the median 
and the dominant contributions in the cross-section. We first determine the scale M^^ that 
gives 50% of the cross- sect ion, and then determine the characteristic transverse scale that 
characterizes the "tail", 80% of the total cross-section ,it's upper boundary will be denoted 
as Mj|. The value of 80% is chosen to quantify that is relevant for the process. The 
more detailed analysis will be given elsewhere. 

Another important scale is the scale of the invariant mass of the dipole pair when 
the partial width for dipole-nucleon scattering reaches one for the central impact parameter 
6 = 0. This scale characterizes the onset of the black disc limit, and is close to that obtained 



in the similar analysis in refs. [18|, ll9|, |20 



As an application of the above discussed results we find that at sufficiently high energies 
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the variety of coherence lengths substitutes the loffe length -l/2mjva; [ll|, 12] familiar from 



QED [13|, [IJl . The derived pattern is different also from the one expected in the application 
of the DGLAP formalism for the scattering off the photon, where the logs come from the 
region of integration over the transverse momenta -C fc^ -C Q^- We pay special attention 
to the coherence length that corresponds to the configurations responsible for the onset of 
the black disk limit (and dominant in the total cross-section). We show that taking into 
account the contribution of dipoles with large masses and of black disc limit to explain the 
probability conservation leads to the coherence length Ic increasing more slowly with the 
energy. This is the because the coherence length for a given process is 

l, = s/{M\s) + Q'), (1.1) 

where M^(s) is the typical M important in the inelastic (total) cross section. Since effective 
is increasing with the energy -see above discussion - the coherence length for the onset 
of the black limit is Ic < (s/(M^(s) + Q^)), and calculations show that this coherence length 
increases with the energy much slower than loffe length, namely as Ic ~ s^'^. 

The increase with energy of the scale has been discussed in detail in js], |l8|, 19, 21 



However, the slowing of the increase of the coherence lengths with energy in the kinematics 
corresponding to the onset of the BD regime, the modified Gribov formulae for the structure 
functions in the BD regime are our new results. 



II. THE DIPOLE APPROXIMATION APPROACH. 

We use in our calculations the dipole approximation, combined with the kt factorization. 
It follows from the QCD factorization theorem 2^ that in the LO approximation inelastic 



cross section of the scattering of the longitudinally polarized virtual photon off a hadron 
target is calculable in terms of the light-cone wave functions of the virtual photon: 
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(2.1) 

Here is the transverse momentum of the constituent within the dipole and kt is the 
transverse momentum of the exchanged gluon. The tensor T"* is the sum of the diagrams 
describing imaginary part of the amplitude for the gluon scattering off the target T. The 
integral of T"'' over (Pkt is proportional to the gluon structure function of the target T. The 
function if) is the wave functions of the virtual photon , and is the photon momentum. 
Within the LO accuracy which we use in the analysis this formula can be rewritten as 



cr 



f^dz J d\{Vim,z)fx'G{x',4r^). (2.2) 



Here the derivatives are over r^t, M = r^/(z{l — z)) is the invariant mass of the dipole, 
x' = (M^ + Q^)/s. The function G(a;',4r^) is the integrated gluon function. The above 
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equation is the generalization of LO DGLAP, and BFKL approximations to the interaction 
of the high energety dipole with the wave function given by pQCD, which accounts for the 
kt factorization theorem. (Within the DGLAP approximation the structure function xG in 
the above formulae depends on Q^,x = Q^/s.) 

The above equation has rather general justification at small x, where the LT approxi- 
mation of pQCD is applicable. Actually it has been understood in preQCD period that in 
a quantum field theory, where the coherence length significantly exceeds the radius of the 
target T at large energies, the electroproduction amplitude in the target rest frame is given 
by a dispersion integral over [26,] . The pQCD guarantees another general property: in the 
target rest frame description: the smaller size of the configuration in the wave function of 
projectile photon leads to a smaller interaction with the target. This property is ensured by 
the kt factorization theorem for the interaction of sufficiently energetic dipole with a target. 
Both properties together give the LO eq l2.2l . In the NLO approximation the structure of 
formulae is the same except the appearance of the additional qqg, ... components in the wave 
function of photon due to the necessity to account for the QCD evolution of the photon wave 
function |6|]. 

For the cross section of the interaction of the longitudinal photon, explicitly differentiating 
the photon wave function, we find: 



Here = 4/3 for the colorless dipoles build of color triplet constituents. The masses 
of the constituents of the dipole were neglected since we restrict our consideration to the 
spatially small dipoles only. 

Similarly, for the pQCD contribution to the cross-section of the transverse photon we 



Here in the practical use we introduce a cut off M^z(l — z) > uf, where Ut is a lower cut 
off, beyond which we can not use the perturbation theory, and in order to obtain the real 
full cross-section we must add the contribution of the aligned jet model (AJM). 

We shall use this equation for the illustrative calculations only, and check it's usefulness by 
studying the cut-off dependence. The key difference from the longitudinal photon case is that 
the asymmetry in the z configurations (the aligned jet model) gives significant contribution 
into the transverse cross-section even at relatively high energies ( this contribution dominates 
in diffractive processes at HERA energies). These contributions correspond to large invariant 
masses even when the actual transverse momenta are small, since = r^/{z{l — z)), 
and large for z —>■ 0, 1. We expect however that for sufficiently large energies the relative 
contribution of the alligned jet model decreases since the contribution of the symmetric 
pQCD configurations increases with energy much more rapidly. Indeed, the contribution of 
the low Tt is multiplied by a structure function at the second argument of order Ar^. 

The formulae discussed above in the momentum space obtain the transparent form in the 
coordinate space description. In the leading order pQCD, the cross-section of the interaction 
of the longitudinal and transverse photons off the nucleon have the form : 




(2.3) 



find: 




(2.4) 




(2.5) 
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Here and below d is the transverse diameter of the dipole, and ijj is the hght-cone wave 
function of the photon and a is the cross section for the quark dipole scattering off target 
T. 

The connection between (P and transverse momenta and invariant masses is given by 
d"^ = 7r^/(4r^), = r?/z(l — z) j23|. The cross-section of the dipole -hadron interaction is 



given by the equation j2J, \2l 

a{d, x') = F\'K^/?,)d^a,{M^)x'G{x', M^). (2.6) 

Here d is the transverse scale of the dipole, and once again x' = {M^ + Q'^)/s. The derivation 
of this formula heavily uses validity of the LO kt factorization theorem -see the above 
discussion. 

The wave functions of the longitudinal and transverse dipole are given by 



\M^,d)\' = {Q/7,^)ae.^.Y.elQ'z\l- zfKl{eh), (2.7) 

1 

where = Q'^z{\ — z) ,and 

n f 

\M^,d)\' = (3/27r2)«e.™.Ee^(^' + (1 - zr){e'Kf{ed). (2.8) 

1 

Here Kq, Ki are the standard zero and first order MacDonald functions. 

Note that with the LO accuracy the cross-section representations as integral over trans- 
verse momenta and over dipole sizes are just connected by a direct and inverse Fourie 
transformation. 

In the practical application of the dipole model for the longitudinal photon the above 
formulae can be further simplified since it may be shown that the dominant contribution 
in the pQCD comes from z=l/2. In this case it is possible to further improve the model, 
introducing a parameter A, such that = A/M^, and determining it from the best fits with 
the data on HERA structure functions and J/ip production. For the symmetric configura- 
tions, taking into account the above connection between d'^ and , we have d'^ = ttVM^. 
However it is possible to show that the results do not change for the whole region of A > 4 



18l | , and thus are insensitive to the precise value of the coefficient in the formula connecting 
d"^ and M^. At high energies a new QCD phenomena appear- the increase with energy of 
the transverse momenta of constituents within the dipole. This property is absent in the 
LO, NLO prescriptions for the calculation of Feynman diagrams of pQCD where the con- 
tribution of the kinematical region: 3> is neglected. In other words, the LO/NLO 
are internally inconsistent at sufficiently large energies since they ignore the basic property 
of pQCD-smaller size of dipole-faster increase of a cross section with energy. We explained 
above that this generalization is necessary to account for basic property of pQCD that the 
interaction of color neutral dipole of zero transverse size with a target should be 0. 

The next practical question is how to parameterize the structure functions for the realistic 
calculations . Here we shall use CTEQ5L distributions |27| . The CTEQ5L distributions has 
been shown to be in good agreement with HERA data [l9[], in the range of x ~ 10~^, 10~^. 
To extrapolate to very small x, we shall use the approximate formulae in the form: 

xG ~ a(M2)/x=(^^'). (2.9) 
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Here the functions 



a{M^) = 2.00123 - 1.69772 ■ 10/M^ + 3.07651/y'MVlO. - 0.228087 • logMVlO., (2.10) 

c(M2) = 0.045 log(M2) + 0.17, (2.11) 

where is in GeV^. This formula is the fit to the observed behavior of the structure 
functions in HERA for 150GeV^ > > 3 GeV^ , made by ZEUS and HI collaborations 
28( 1 . We shall use this function also in the kinematics where the CTEQ parton distributions 
are absent. Note that observed increase of c(M^) with can be explained as due to 



evolution of parton distributions cf. ref.[2^ and within the resumation approach [ol, 



III. THE HARD FRAGMENTATION PROCESSES. 

We have carried both analytical and numerical analysis of the transverse scales. The 
analytical calculation has been made for the toy model based on the double logarithmic 
approximation and will be published in the full version of the paper. The numerical results 
will be presented below. We shall extrapolate our results to energies of order s ~ 10^ 
GeV^. These energies are unattainable due to the fact that there is no e-p collider for such 
energies. (The proposed e-p DIS facility at the LHC may reach the energies of order 10^ 
GeV^. However these results are interesting both from the theoretical point of view (probing 
the limits of the pQCD at very high energy), and from the practical point of view, giving 
the information about the parton distributions at the LHC. The relation of our results to 
the processes at the LHC will be discussed in section V. 



A. Numerical analysis: longitudinal photons. 

Let us analyze the integrand in the integral representation of the cross-section 12.51 as 
the integral over (P. This integral can be rewritten as integral over = 4:/rf, = 
{7i'^/4:(P) jl^l, and we use the dominance of the symmetric configurations z = 1/2. This 
approximation is good since z{l — z) is a slow function of z for a large range of z. With 
the logarithmic accuracy such integral must numerically be the same as the integral over 
in the momentum representation of the corresponding cross-section. We carried our 
calculations both in the coordinate and momentum representations and found that the 
results coincide within a 10% accuracy. 

Let us first consider the quantity da^/dM'^. This quantity in the z=l/2 approximation 
is proportional to the transverse momentum rt jet distribution in the fragmentation region. 

In order to characterize the behavior of the density we calculate two scales: Mf and M^. 
The first of these scales characterizes the maximum of the density. This scale lies Mf < Q^, 
for all the values of the energy less than 10^^ GeV^. In line with the expectations based 
on the toy model we see that this maximum however slowly increases with energy for every 
value of Q^, with the energy dependence for > 5 GeV^ being approximately s*^"^. 

The scale Mf however does not give a full characterization of the density, the reason is 
that the rapid decrease of the square of the derivative of the wave function is compensated 
partly by the rapid increase of the structure function with the increasing M^. In order to 
characterize this effect it is natural to determine the new scale that corresponds to the 
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cut off in that gives, say, 50% of the total cross- sect ion, and another scale M^2 that 
corresponds to 80% of the cross-section and characterizes the magnitude of the tail. 

The determination of M^^ is given in Table 1 for the characteristic virtualities = 
5,10,20,40,100 GeV^ and for reahstic energies (up to those achievable at LHC). We see 
that for not very small (starting from ~ 10 GeV^) and for HERA energies this scale 
lies beyond 0.75Q^, where is an external virtuality, thus justifying the use of DGLAP at 
these energies. We see from the tables that this median scale rapidly increases with energy, 
and will overcome at = 10 — 20 GeV^ already at LHC energies and may significantly 
overcome if we shall go beyond LHC. 

Let us consider the second scale, Mj|. We see that the tail exists even at HERA energies 
for < 200 GeV^, when the M^l exceeds Q^. The tail also increases with the increase of 
energy, even more rapidly then the median scale M^^. 

The scale is always larger than Q^. If we extrapolate to the very large (though 
unrealistic) energies we see that it exceeds for very small ~ 5GeV^ by a factor of 
5. In other words, for sufficiently small momenta and high energies the leading logarithmic 
approximation leads to the very long tails, that are strictly beyond the control of usual 
prescriptions for DGLAP and BFKL approximations and must be taken into account using 
kt factorization theorem, at least in the LO approximation. 

We conclude that for sufficiently small external virtualities and high energies (large 1/x) 
one cannot use the naive DGLAP approximation where transverse momenta in the upper 
rung of the ladder are -C Q, while the the dipole approximation provides a reasonable 
description of the cross- sect ions. In this approximation the LO DGLAP/resummed ladder 
gives the cross-section of the dipole scattering of the target. 

B. The evaluation of the transverse scale: transverse photons. 

The same two scales and Mf that we defined for longitudinal photon can be defined 
for the transverse photon. 

In this case however the use of the pQCD expression poses a problem due to a poten- 
tial large contribution of the nonperturbative physics, described by the aligned jet model. 
Another difference from the longitudinal photon case is the large contribution of the asym- 
metric configurations, with z 1/2. These contributions are dominant in the AJM model, 
and lead to large invariant masses even for small transverse momenta, due to relation 
= rj/{z{l - z). In this talk we are concerned only with the contribution of the pQCD. 
In order to exclude the contribution of the AJM (nonperturbative QCD) we impose a cutoff 
M'^z{l — z) > A^, which removes the contribution of the low transverse momenta. The use of 
the pQCD formulae, even for illustrative purposes, would be justified only if the cross section 
depends weakly on the cutoff. We see however, that for < 20 GeV^ even for the energies 
of the order 10^ GeV^ there is a strong dependence (of the order 20%) of the maximum of 
the curve da/dM'^ on a very small change of the cutoff (from A = 0.75 to A = 1 GeV^. The 
data for the median and the tail transverse scales M^^, M^2 are presented in Tables 3, 4. We 
see that for HERA energies there is a dependence of order 10% of these scales on the cut 
off even at ~ 100 GeV^. This dependence however decreases with energy. 

The detailed analysis of the relative contribution of the pQCD and aligned jet model 
will be given in ref. j3l|. In tables 4, 5 we presented our results for transverse scales at 
Q'^ = 20,40, 100 GeV^. We considered the integrand for cross-section as a function of the 
invariant mass and cut off, after integrating over permitted z for given invariant mass and 
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cut off. 

The scale Mf < for all energies and > 10 GeV^ behaves similar to the longitudinal 
photon, However this scale is shifted to smaller invariant masses relative to longitudinal 
photons. For every given virtuality the scale continues to increase with the increase of 
energy, with approximately the same rate as for longitudinal photons. 

The characteristic feature of the invariant mass distribution for the transverse photons 
is however that it becomes much broader, and increases with energy more rapidly than for 
longitudinal photons. As a result although the transverse scales for HERA energies are 
slightly smaller than for longitudinal photons, for high energies, achievable at LHC they 
already overcome them. The rate of increase is ~ s"'^ for LHC energies and ~ 50 
GeV2. 

It is once again instructive to look where = as a function of energy. We see that 
the tail starts (with logarithmic accuracy) for higher than for the case of the longitudinal 
photons. 

Since our results show qualitative stability with the change of the cut off and the depen- 
dence of cut off (and AJM contribution ) decrease with and energy, we expect that our 
results remain qualitatively the same if we shall take into account the AJM contribution 
explicitly. 

Let us note that the appearance of the tails at sufficiently large energies (small x) can 
be expected from the comparison with the Gribov formulae for the BD regime [26']. Indeed, 
assuming a smooth matching of the pQCD and BD regime the perturbative distribution 
must match the black limit spectral density , and this density leads to the tail with masses 
increasing with energy. Thus the existence of a large mass tail beyond the naive DGLAP 
approximation seems to be a necessary property of pQCD near the black limit. 

Finally, we expect that our results will not change in the NLO/resumed model. The 
reason is that NLO effects are partly taken into account in the gluon distribution fits to 
experimental data. 

C. The evaluation of the black scale. 

We can now combine the constraint for the partial wave at zero impact parameter to 
become unit with the CTEQ5 structure function to determine the invariant masses that 
correspond to the onset of the black limit. 

The typical dependence of the black limit onset scale on energy is presented in Table 
3. for the gluonic dipole. We consider two cases: one when the partial wave F at the central 
impact parameter reaches 1, another when it reaches 1/2 [3|. Indeed, when the partial wave 
reaches 1/2 the probability of inelastic interactions reaches 3/4, i.e. interactions become 
strong and pQCD can not be used any more [3|. 

The characteristic transverse momenta that corresponds to the black scale for the gluonic 
dipole for the energies s = 10^ — 10^ GeV^ are 2 — 2.5 GeV/c if we impose F = 1 condition 
or 4 — 4.5 GeV/c if we impose F = 1/2. Note that for HERA energies (s = 10'^ GeV^ we do 
not obtain the black limit for gluonic dipole for > 40 GeV^, while for lower values of 
we obtain formally A;* ~ 1 GeV, the value that seems beyond the limits of the applicability 
of the method we used here to determine this scale. 

For fermionic dipoles we see that there is no black limit for HERA energies, although 
nonperturbative effects (corresponding to F = 1/2) seem to start to appear for kt < 1 GeV. 
For the energies s ~ 10® — 10^ GeV^ the transverse scale is ~ 2 — 2.5 GeV for the partial 
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wave r = 1 condition, 3 — 3.5 GeV respectively for F = 1/2 conditions. These results are in 
good agreement with the previous determination of these scales 0, 18, 20 . 



For our purposes it will be important to determine the rate of increase of the M"^. This 
rate does not depend on the type of the dipole (fermionic or gluonic) or on the partial wave 
condition. 

For the gluonic dipole in DIS for s up to s ~ 10^ GeV^ the black limit scales increase 
relatively slow as s°'^ . However at the energies above s ~ 10'' GeV^ the increase speeds up 
with ~ s^''^ starting from s ~ 10* GeV^ energies ( for the exponent for partial wave 1/2 
we obtain 0.39, for the exponent with partial wave equal to 1 we obtain 0.38), and small 

< 40 GeV^. Note that this region is where we expect that there is no dip influence. 

We show that for moderately large Q^, where the coupling constant is small the cross- 
section will be dominated by the black limit. 

We found from this subsection that for the DIS processes at very high energies there exist 
3 regimes. First the nonperturbative black disk regime. This regime is valid for invariant 
masses up to M^, then the standard pQCD regime described by DGLAP and then the new 
pQCD regime, where one can not use DGLAP, BFKL resummed models directly, but must 
combine them with the dipole approximation. In this new regime a virtual photon creates 
a parton pair with invariant mass larger than the external virtuality, that scatters off the 
target by the DGLAP/resummed ladder. Thus in the s— plane we find three regions: the 
black limit region-the virtualities where the cross-section is dominated by the black limit, 
the region where the cross-section is dominated by the usual pQCD, and the new regime 
area where the tail gives a significant contribution to cross-section. 



D. A hadron structure function in the black disc limit 

The above formulae make possible first quantitative evaluation of the numerical coefficient 
in the asymptotical expression for the structure function of DIS in the black disc limit 
suggested in ref . jl, 0] . 

The asymptotic behavior of the hadron structure functions at s — oo in QCD is rather 
close to Gribov formulae for the total cross section of DIS derived long ago for heavy nuclear 
target in the black disc limit 2^. In the case of the transverse photon 



Fs = 27tR^{s)Q'^kp{M^ oo)(ln(a;o/x)), (3.1) 
where p is the spectral density for the transition e'^e~ hadrons : 

p(M^) = "^f f ^ ^"^^""^^ (3.2) 

The density p{M'^) tends to constant at large invariant masses. Rt is the maximal impact 
parameter, where the partial wave reaches 1. Fourie transform of two gluon form factor 
of a nucleon: F = 1/(1 + / p'^Y measured in the hard diffractive processes at HERA 
and at FNAL describes dependence of a partial amplitude / on the impact parameter b : 
f{b,x) oc {1/ x)'^exp~^^ . Here p^ ~ IGeV^ . In the LT approximation the structure function 
of a hadron in DIS is usually parametrized as F2 ~ (1/a;)^. Finally we obtain: 

Rt = \{Ml)ln{xQ/x)/p (3.3) 
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Here xq ~ 10~^, A slightly more complicated reasoning which uses completeness over the 
hadronic states shows that Eq J3.1l is valid even if to account for the gluons and quark- 
antiquarks pairs in the wave function of the dipole. The coefficient n in eq. I3.1l is determined 
by the upper limit of the integration over and is equal to ~ s'^. Our results show 
that this coefficient weakly depends on energy, with k, = 0.3 — 0.4 at the LHC energies and 
beyond. 



IV. THE COHERENCE LENGTH 

In the previous sections we determined the absolute values and the energy dependencies 
of the effective transverse scale and the black limit scale at high energies which allowed us 
to evaluate coherence length. The coherence length corresponds to the time, such that the 
dipole fluctuation exists this time interval at a given energy. The original suggestion of the 
existence of the coherence length in the deep inelastic scattering was first made by loffe, 
Gribov and Pomeranchuk 11, 2^ It was found already in the sixties by loffe 12 1 that the 



coherence length at moderate is / ~ l/2mNXB i-e. it linearly increases with energies. 
At higher energies we have a whole variety of local coherence lengths. Here we consider the 
coherence length, that corresponds to configurations near the onset of the black limit that 
dominate the cross-section at very high energies. 

Since near the onset of the black limit ~ we conclude, that in the vicinity of 
black limit i.e. for sufficiently high energies this coherent length increases like Ic ~ s^'^, i.e. 
it increases with energy much slower than for lower energies, when loffe result is applicable. 

The important characteristic of the hard processes is the local gluon density produced 
in the hard high energy process. Rough estimate of the three dimensional space density of 
gluons (number of gluons in a unit of the space volume) gives 

n ~ N/ikn/r',) > 4r,2)/(U/r2), (4.1) 

where N is a total number of the emitted gluons. 

At the LHC energies the coherence length increases with energy as ~ s^'^ for sufficiently 
small (~ s^'^^ near the energies of order 10^^ GeV^ ). Using known dependence of the 
gluon structure function on energy, we obtain that the number of gluons increases like 
while the same formulae extrapolated to superhigh energies gives This means, that 
the naive estimate gives slowly decreasing gluon number density n ~ s"'''^ — s~^'^. Here, 
however, we neglected the increase of the transverse momenta that may also contribute to the 
increase of a local gluon density. The realistic estimate of this increase is however impossible 
in the LO approximation. Consequently, within our accuracy we cannot determine whether 
the local (three-dimensional) gluon density is increasing or slowly decreasing. To address 
this question more quantitatively we have to go beyond the naive estimate and perform a 
detailed calculation of the local density, taking into the account the inhomogeneity of the 
ladder. Such analysis will be done elsewhere. 

We conclude that there are indications that the local gluon density may appear large in 
the black disc limit regime. 

The opposite example, of high light cone densities and low local density rapidly decreasing 
with energy was considered some time ago by Mueller in QED [30]. However in that case 
there is no rapid increase of interaction with energy. 

It is worth noting that discussed above pattern for the energy dependence of the coherence 
length leads to a change of the structure of the fast hadron wave function as compared to the 
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Gribov picture where the longitudinal size of the hadron is determined by the wee parton 
cloud and energy independent ~ where /i is the soft scale. On the other hand a slower 
rate of the increase of the coherent length with energy than I/ttinx leads to decrease of 
the longitudinal size of the hadron with energy. The typical size is determined by the BD 
momentum at a given impact parameter for a particular energy. Moreover since the BD 
momentum is larger for small impact parameters the nucleon has a form of a concave lens. 
It is of interest also that for the zero impact parameter the longitudinal size of a heavy 
nucleus is smaller than for a nucleon. 



V. EXPERIMENTAL CONSEQUENCES 

The current calculations of cross-sections of hard processes at the LHC are based on the 
use of the DGLAP parton distributions and the application of the factorization theorem. 
Our results imply that the further analysis is needed to define the kinematic regions where 
one can use DGLAP distributions. We showed in the paper that for DIS at high energies 
there are kinematic regions where one is forced to use a kt factorization and the dipole model 
instead of direct use of DGLAP. A similar analysis must be made for the pp collisions at 
LHC. This analysis is however more complicated since the proton can not be approximated 
by a dipole and thus the DIS results for the same energies and external virtualities can not 
be transferred directly to pp case. This is a problem for a future work. 

The hard processes initiated by the real photon can be directly observed in the ultraphe- 
ripheral collisions [20|. The processes where a real photon scatters on a target, and creates 
two jets with an invariant mass M^, can be analyzed in the dipole model by formally putting 

= 0, while is an invariant mass of the jets. Then with the good accuracy the spectral 
density discussed above will give the spectrum of jets in the fragmentation region. Our 
results show that the the jet distribution over the transverse momenta will be broad with 
the maximum moving towards larger transverse momenta with increase of the energy and 
centrality of the 7^4 collision. 

Finally, our results can be checked directly, if and when the LHeC facility will be built 
at CERN. 
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Table 1. The scale M^^ (50% of the total cross-section) for longitudinal photons in DIS 





Q2 = 5 GeV^ 


= 10 GeV2 


= 20 GeV2 


— 40 GeV^ 


— 100 GpV^ 


s = W GeV2 


6.5 GeV^ 


9 GcV^ 


16 GeV2 


28 GpV^ 


60 GpV^ 


s = 10-' GeV^ 


7 GcV2 


10.5 GcV^ 


17.5 GeV^ 


31 GpV^ 

O ± VJC V 


68 GpV^ 


s = 106 GeV2 


7.5 GeV^ 


11 GeV^ 


19 GeV2 


34GpV2 


77GpV2 


s = 10^ GeV2 


8 GeV^ 


12 GeV2 


21 GeV2 


37 GeV^ 


87 GeV^ 


Table 2. The scale (8C 


)% of the cross-section ] 


for lonEfituf 


inal nhotons 




g2 ^ 5 GeV^ 


= 10 GeV^ 


= 20 GeV^ 


= 40 GeV^ 


= 100 GeV^ 


s = 10^ GeV^ 


10.5 GcV2 


17 GeV2 


30 GeV2 


54 GeV2 


120 GeV^ 


s = 10^ GeV^ 


11.5 GcV2 


19 GeV2 


34 GeV2 


60 GeV^ 


140 GeV2 


s = 10^ GeV^ 


12.5 GeV^ 


21 GeV2 


38 GeV2 


67GeV2 


160GeV2 


s = 10^ GeV2 


14 GeV^ 


23 GeV2 


42 GeV2 


75 GeV^ 


180 GeV2 


Table 3. The scale M^^ i 


or transverse photons in DIS. 



s = 10^ GeV2 



s = 10^ GeV^ 



s = 10^ GeV^ 



s = 10^ GeV2 



g2 = 20 GeV2 



24(20) GeV^ 



27(23) GeV^ 



32(28) GeV^ 



38(33) GeV^ 



40 GeV2 



31(26) GeV^ 



37(33) GeV^ 



43(38) GeV^ 



52 (48) GeV^ 



Table 4. The scale M^2 for transverse photons in DIS 



100 GcV^ 



50(45) GeV2 



59 (55)GeV2 



73 (70)GeV^ 



90 (86) GeV^ 





g2 = 20 GeV2 


= 40 GeV^ 


= 100 GeV^ 


s = 10^ GcV2 


73 (65) GeV2 


100 (90) GcV2 


160 (150) GcV^ 


s = 10^ GeV^ 


83 (75) GeV2 


120 (110) GeV^ 


200 (190) GeV^ 


s = 10^ GeV2 


96 (88) GeV^ 


140 (130) GeV2 


260 (250) GeV2 


s = 10^ GeV^ 


110 (100) GeV2 


180 (170) GeV2 


330 (320) GeV2 



Table 5. The scale kt for the onset of the black dis^ 





r = 1 


r = 1/2 


s = W GcV2 


1 GcV 


1.6 GcV 


s = 10^ GcV2 


1.6 GeV 


2.3 GeV 


s = 10^ GcV^ 


2.3 GeV 


3.2 GeV 


s = W GeV^ 


3.3 GeV 


4.5 GeV 



c regime. 



The values of here in the tables 3,4 correspond to cut offs 1 GeV, and 0.75 GeV 
(in the brackets) 
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